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Abstract. In the oblivious buy-at-bulk network design problem in a 
graph, the task is to compute a fixed set of paths for every pair of source- 
destinations in the graph, such that any set of demands can be routed 
along these paths. The demands could be aggregated at intermediate 
edges where the fusion-cost is specified by a canonical (non-negative con- 
cave) function /. We give a novel algorithm for planar graphs which is 
oblivious with respect to the demands, and is also oblivious with respect 
to the fusion function /. The algorithm is deterministic and computes the 
fixed set of paths in polynomial time, and guarantees a 0(log n) approx- 
imation ratio for any set of demands and any canonical fusion function 
/, where n is the number of nodes. The algorithm is asymptotically op- 
timal, since it is known that this problem cannot be approximated with 
better than fi(logn) ratio. To our knowledge, this is the first tight anal- 
ysis for planar graphs, and improves the approximation ratio by a factor 
of log n with respect to previously known results. 



1 Introduction 

A typical client-server model has many clients and multiple servers where a 
subset of the client set wishes to route a certain amount of data to a subset of 
the servers at any given time. The set of clients and the servers are assumed 
to be geographically far apart. To enable communication among them, there 
needs to be a network of cables deployed. Moreover, the deployment of network 
cables has to be of minimum cost that also minimizes the communication cost 
among the various network components. This is what we roughly call as a typical 
network design problem. The same problem can be easily applied to many similar 
practical scenarios such as oil/gas pipelines and the Internet. 

The "Buy-at-Bulk" network design considers the economies of scale into ac- 
count. As observed in [5J, in a telecommunication network, bandwidth on a 
link can be purchased in some discrete units u\ < U2 < ■ ■ ■ < Uj with costs 
ci < C2 < ■ ■ ■ < Cj respectively. The economies of scale exhibits the property 
where the cost per bandwidth decreases as the number of units purchased in- 
creases: c\ju\ > C2/112 > ...Cj/uj. This property is the reason why network 
capacity is bought/sold in "wholesale", or why vendors provide "volume dis- 
count" . 
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The generalized form of the buy-at-bulk problem is where there are multiple 
demands from sources to destinations, and it is commonly referred as Multi-Sink 
Buy-at-Bulk (MSBB). Typically, the demand flows are in discrete units and are 
unsplittable (indivisible), i.e., the flow follows a single path from the demand 
node to its destination. These problems are often called "discrete cost network 
optimization" in operations research. 

As mentioned in 16. , if information flows from x different sources over a 
link, then, the cost of total information that is transmitted over that link is 
proportional to f(x), where / : Z + — > M + . The function / is called a canonical 
fusion function if it is concave, non-decreasing, /(0) = and has the subadditive 
property f(xi +X2) < f{x\) + f{x2), Wxi f X2, {x\ + x 2 ) £ Z + . Generally, MSBB 
problems use the subadditive property to ensure that the 'size' of the aggregated 
data is smaller than the sum of the sizes of individual data. 

We study the oblivious buy-at-bulk network design problem (MSBB) with 
the following constraints: an unknown set of demands and an unknown concave 
fusion cost function /. An abstraction of this problem can be found in many 
applications, one of which is data fusion in wireless sensor networks where data 
from sensors is aggregated over time in multiple sinks. Other application include 
Transportation & Logistics (railroad, water, oil, gas pipeline construction) etc. 
Many of these problems are formulated as networks on a plane that can be 
mapped to planar graphs. 

1.1 Problem Statement 

Assume that we are given a weighted graph G = (V,E,to), with edge weights 
to : E — > Z+. We denote ro e to be the weight of edge e. Let di — (si,ti) be a 
unit of demand that induces an unsplittable unit of flow from source node Si € V 
to destination node ti £ V. Let A = {d\,d 2) ... ,d r } be a set of demands that 
are routed through paths in G. It is possible that some paths may overlap. The 
flow of these demands forms a set of paths P(A) = {p(di),p(d,2), . . . ,p(d r )}. 

There is an arbitrary canonical function / at every edge where data aggre- 
gates. This / is same for all the edges in G. Let f e (A) — {p{di) : e £ p(di)} 
denote the set of paths that use an edge e £ E. Then, we define the cost of 
an edge e to be C e (A) — f(\<p e (A)\) ■ w e . The total cost of the set of paths is 
defined to be C(A) = J2 e eE^e(A). For this set A, there is an optimal set of 
paths P* (A) with respective cost C* (A) . The approximation ratio for the paths 
P(A) is defined as ' M£>E$E5 optimization problem on input A is to 

find a set of paths P(A) that minimizes the approximation ratio. We note that 
MSBB is NP-Hard as the Steiner tree problem is its special case (when f(x) = 1 
and when there is only one destination node) [151 . 

An oblivious algorithm A u for the MSBB problem, computes a fixed set of 
paths, denoted P(G) for every pair of source destination nodes in V. Given any 
set of demands A, the path p(di) for each di = (si,U) £ A, is the fixed path in 
P(G) from Si to ti. This gives a set of paths P(A) to route the demands A. We 
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define the approximation ratio of A u, as: 

C(A) 



A.R.(Aobi) — max 



a C*(A)' 

We aim to find algorithms that minimizes the above approximation ratio for 
any canonical function / which is unknown to the algorithm. The best known 
oblivious algorithm is by by Gupta et al. [H] and provides approximation ratio 
0(log 2 n) for general graphs. No better result is known for planar graphs. This 
problem is NP-hard, since MSBB is NP-hard. 



1.2 Contribution 

We provide an oblivious algorithm Find Paths for MSBB problems in planar 
graphs. Our algorithm is deterministic and computes in polynomial time a fixed 
set of paths that guarantees 0(log n)-approximation ratio for any canonical func- 
tion / (where / is unknown to the algorithm) . We also give a lower bound for the 
approximation ratio for A to be of J?(logn), where n is the number of nodes in 
the graph. A lower bound of X2(log n) for planar graphs is provided in the context 
of the online Steiner tree problem by Imase and Waxman [10] . Thus, our bound 
is tight with respect to planar graphs. It is also a logn factor improvement over 
the best previously known result [9]. 

We build the set of paths based on sparse covers (see [14] for an overview 
of sparse covers). A 7-cover consists of clusters where for each node there is 
some cluster that contains its 7-neighborhood. We construct O(logn) levels of 
covers with exponentially increasing locality parameter 7. For every cluster we 
elect a leader. For any pair of nodes u, v we identify an appropriate common 
lowest-level cluster that contains both u and v, and the cluster has a respective 
common leader £. Then the path from u to v is formed by connecting successive 
path segments emanating from both u and v and using intermediate leaders of 
lower level clusters until the common leader t is reached. 

In the analysis, we introduce the notion of coloring sparse covers, where two 
clusters that are close receive different color. We show the existence of a sparse 
cover with constant coloring (based on the sparse covers in [4]). This enables us 
to obtain optimal approximation at every level. When we combine all the levels, 
we get an O(logn) approximation. 

1.3 Related Work 

Oblivious Network Design Below, we present the related work on oblivious 
network design and Table [T] summarizes some results and compares our work 
with their's. What distinguishes our work with the others' is the fact that we 
provide a set of paths for the MSBB problem while others provide an overlay 
tree for SSBB version. 

Goel et al. in [6 build an overlay tree on a graph that satisfies triangle- 
inequality. Their technique is based on maximum matching algorithm that guar- 
antees (1+log fc)-approximation, where k is the number of sources. Their solution 
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Tabic 1: Our results and comparison with previous results for data- fusion 
schemes, n is the total number of nodes in the topology, k is the total num- 
ber of source nodes. 



Related 


Algorithm 


Graph 


Oblivious 


Oblivious Approx 


Work 


Type 


Type 


Function / Sources 


Factor 


Lujun Jia 
et al. [12] 


Deterministic 


Random 
Deployment 


X 


/ 


O(logn) 


Lujun Jia 


Deterministic 


Arbitrary Met- 
ric 


X 


/ 


0( log4 ™ ) 

\ log log(n) / 


et al. rn] 


Deterministic 


Doubling Met- 
ric 


X 


/ 


0(log(n)) 


Ashish Goel 
et al. [5] 


Randomized 


General Graph 
A-inequality 


/ 


X 


O(logfc) 


Ashish Goel 
et al. [7] 


Probabilistic 


General Graph 


/ 


X 


O(l) 


Anupam Gupta 


Randomized 


General Graph 


/ 


/ 


OQog'n) 


et al. [9] 


Randomized 


Low Doubling 


/ 


/ 


O(logn) 


This paper 


Deterministic 


Planar 


/ 


/ 


O(logn) 



is oblivious with respect to the fusion cost function /. In a related paper [7], Goel 
et al. construct (in polynomial time) a set of overlay trees from a given general 
graph such that the expected cost of a tree for any / is within an 0(l)-factor of 
the optimum cost for that /. 

Jia et al. in [T2] build a Group Independent Spanning Tree Algorithm (GIST) 
that constructs an overlay tree for randomly deployed nodes in an Euclidean 2 
dimensional plane. The tree (that is oblivious to the number of data sources) 
simultaneously achieves 0(logn)-approximate fusion cost and 0(l)-approximate 
delay. However, their solution assumes a constant fusion cost function. We sum- 
marize and compare the related work in Table [TJ 

Lujun Jia et al. provide approximation algorithms for TSP, Steiner Tree 
and set cover problems. They present a polynomial-time (0(log(n)), 0(log(n)))- 
partition scheme for general metric spaces. An improved partition scheme for 
doubling metric spaces is also presented that incorporates constant dimensional 
Euclidean spaces and growth-restricted metric spaces. The authors present a 
polynomial-time algorithm for Universal Steiner Tree (UST) that achieves poly- 
logarithmic stretch with an approximation guarantee of 0(log 4 nj loglog(n)) for 
arbitrary metrics and derive a logarithmic stretch, 0(log(n)) for any doubling, 
Euclidean, or growth-restricted metric space over n vertices. They provide a 
lower bound of f2(logn/loglogn) for UST that holds even when all the vertices 
are on a plane. 

Gupta et al. [5] develop a framework to model oblivious network design prob- 
lems (MSBB) and give algorithms with poly-logarithmic approximation ratio. 
They develop oblivious algorithms that approximately minimize the total cost 
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of routing with the knowledge of aggregation function, the class of load on each 
edge and nothing else about the state of the network. Their results show that if 
the aggregation function is summation, their algorithm provides a 0(log 2 n) ap- 
proximation ratio and when the aggregation function is max, the approximation 
ratio is 0(log 2 n log log n). The authors claim to provide a deterministic solution 
by derandomizing their approach. But, the complexity of this derandomizing 
process is unclear. 

Non-Oblivious Network Design There has been a lot of research work in 
the area of approximation algorithms for network design. Since network design 
problems have several variants with several constraints, only a partial list has 
been mentioned in the following paragraphs. 

The "single-sink buy-at-bulk" network design (SSBB) problem has a single 
"destination" node where all the demands from other nodes have to be routed 
to. Network design problems have been primarily considered in both Operations 
Research and Computer Science literatures in the context of flows with concave 
costs. The single-sink variant of the problem was first introduced by Salman et al. 
|15j . They presented an 0(log n)-approximation for SSBB in Euclidean graphs by 
applying the method of Mansour and Peleg [13] . Bartal's tree embeddings [2] can 
be used to improve their ratio to O(lognloglogn). A (9(log 2 n)-approximation 
was given by Awerbuch et al. [1] for graphs with general metric spaces. Bartal 
et al. [3] further improved this result to O(logn). Guha [8] provided the first 
constant-factor approximation to the problem, whose ratio was estimated to be 
around 9000 by Talwar [16]. 

Organization In the next section, we present some definitions and notations 
used throughout the rest of the paper. Section [3] describes the concept of sparse 
covers and shortest path clustering. In addition, we derive a coloring for the 
cover. In section 2] we describe Fixed Paths and Find Paths algorithms that build 
a set of shortest paths between all pairs of nodes in G. Section [5] provides the 
analysis of the Find Paths algorithm as well as the main theorem of this paper. 
Finally, we discuss our contribution and future work in section [6] 

2 Definitions 

Consider a weighted graph G — (V,E,tt>), where to : E — > Z + . For any two 
nodes u, v G V, their distance dist(u,u) is the length of the shortest path that 
connects the two nodes in G. We denote by Nk(v) the k -neighborhood of v which 
is the set of nodes distance at most k from v. For any set of nodes S C V, we 
denote by Nk(S) the fc-neighborhood of S which contains all nodes which are 
within distance k from any node in S. 

A set of nodes X C V is called a cluster if the induced subgraph G(X) is 
connected. Let Z = {Xx, X 2 , . . . , Xk} be a set of clusters in G. For every node 
v G G, let Z(v) C Z denote the set of clusters that contain v. The degree of v in 
Z is defined as (3 V (Z) = \Z(v)\, which is the number of clusters that contain v. 
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The degree of Z is defined as f3(Z) = max„gy (3 V (Z), which is largest degree of 
any of its nodes. The radius of Z is defined as rad(Z) = maxxez(^8i,d(X)). 

Consider a locality parameter 7 > 0. A set of clusters Z is said to ^-satisfy 
a node v in G, if there is a cluster X E Z, such that the 7-neighborhood of v, 
Nry(v), (nodes within distance 7 from v) is included in X, that is, N 7 (v) C X. 
A set of clusters Z is said to be a j-cover for G, if every node of G is 7-satisfied 
by Z in G. The stretch o~(Z) of a 7-cover Z is the smallest number such that 
rad(Z) = a(Z) ■ 7. 

We define the following coloring problem in a set of clusters Z. We first 
define the notion of the distance between two clusters X{, Xj G Z, Xi 7^ Xj. 
We say that dist(A i , Xj) < k, if there is a pair of nodes u G X, t and v G Xj 
such that u is fc-satisfied in Xi, v is fc-satisficd in Xj, and dist(u, v) < k. A valid 
distance- /c coloring of Z with a palette of x colors [1,%]) is an assignment of an 
integer color(A) G [l,x] to every X G Z, such that there is no pair of clusters 
Xi, Xj G Z, Xi ^ Xj, with dist(Xi, Xj) < k which receive the same color. The 
objective is to find the smallest x that permits a valid distance-fc coloring. 

3 Sparse Cover 

A 7-cover is sparse if it has small degree and stretch. In [4| Section 5] the au- 
thors present a polynomial time sparse cover construction algorithm Planar- 
Cover(G, 7) for any planar graph G and locality parameter 7, which finds a 
7-cover Z with constant degree, (3 < 18, and constant stretch, a < 24. Here, 
we show that this cover also admits a valid distance-7 coloring with a constant 
number of colors x < 18- 

For any node v G G, we denote by depths (G) the shortest distance between 
v and an external node (in the external face) of G. We also define depth(G) = 
max^gy depth tI (G). The heart of sparse cover algorithm in [H Section 5] concerns 
the case where depth(G) < 7 which is handled in Algorithm Depth-Cover(G, 7). 
The general case, depth(G) > 7, is handled by diving the graph into zones of 
depth 0(7), as we discuss later. So, assume for now that depth(G) < 7. 

The Algorithm Depth-Cover(G, 7), relies on forming clusters along shortest 
paths connecting external nodes (in the external face) of G. For every short- 
est path p, Algorithm Shortest-Path-Cluster(G,p, 47) in [4j Section 3] returns a 
set of clusters around the 47 neighborhood of p with radius at most 87 and 
degree 3. Then, p and all its 27-neighborhood is removed from G producing a 
smaller subgraph G' (with possibly multiple connected components). The algo- 
rithm proceeds recursively on each connected component H of G' by selecting 
an appropriate new shortest path p' between external nodes of H . The algorithm 
terminates when all the nodes have been removed. The initial shortest path that 
starts the algorithm consists of a single external node in G. The resulting 7- 
cover Z consists of the union of all the clusters from all the shortest paths. The 
shortest paths are chosen in such a way that a node participates in the clustering 
process of at most 2 paths, and this bounds the degree of the 7-cover to be at 
most /3 < 6, and stretch s < 8. 
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The analysis in [4j Section 5.1.1] of Algorithm Depth-Cover relies on repre- 
senting the clustering process of G as a tree T as we outline here. Each tree node 
w G T represents a pair w = (G(w),p(w)) where G(w) is a planar subgraph of G 
that is to be clustered, and p(w) is a shortest path between two external nodes 
of G(w). The root of the tree is r — (G,v), where v is a trivial initial path 
with one external node v € G. The children of a node w € T are all the nodes 
w' = (G(w'),p(w')), such that G(w') is a connected component that results from 
removing p(w) and its 27-neighborhood from G. 

Next, we extend [H Lemma 3.1] to show that we can color the clusters ob- 
tained by a shortest path clustering using a constant number of colors. The proof 
of the following three lemmas is given in the appendix. 

Lemma 1. For any graph G, shortest path p G G, the set of clusters returned 
by Algorithm Shortest-Path-Cluster(G,p, 4j) admits a valid distance-"/ coloring 
with 3 colors. 

We can obtain a coloring of Z by coloring the respective levels of the tree T. 
Assume that the root is at level 0. 

Lemma 2. The union of clusters in any level i > of the tree T , admits a valid 
distance-"/ coloring with 3 colors. 

Lemma 3. Algorithm Depth-Cover(G, 7) returns a set of clusters Z which ad- 
mits a valid distance-"/ coloring with 6 colors. 

We are now ready to consider the case depth(G) > 7. Algorithm Planar- 
Cover(G, 7) decomposes the graph into a sequence of bands, such that each band 
Wi has depth 7. The bands are organized into zones, such that zone Si consists 
of three consecutive bands Wj_i , Wi , Wj+i . Thus, zone Si overlaps with bands 
Si-2, S^ 1, Si+i and Si+2- The algorithm invokes Depth-Cover(Si, 37— 1) for each 
zone giving a 7-cover Z with degree j3 < 3 • 6 = 18 and stretch a < 3 • 8 = 24. 

We can obtain the following coloring result. Using Lemma [31 for every zone 
Si we can get a valid distance- (37 — 1) coloring with a palette of 6 colors. This 
implies that we can obtain a valid distance-7 coloring for the zone with at most 
6 colors. Zones Si and Si+3 do not overlap and any two nodes satisfied in them 
(one from each zone) with respect to G must be more than 7 distance apart. 
Therefore, we can color all the zones with three different palettes each consisting 
of 6 colors, so that zone Si, uses the ((« mod 3) + l)th palette. The coloring can 
be found in polynomial time. Therefore, we obtain: 

Theorem 1. Algorithm Planar-Cover(G, 7) produces a set of clusters Z which 
has degree j3 = 18, stretch a < 24, and admits a valid distance-"/ coloring with 
X = 18 colors. 

4 Algorithm 

We describe how to find the shortest path between a given pair of nodes in a 
graph G = (V, E, ro). To find such paths, we use Algorithm Find Paths (Algorithm 
[2]) which relies on Algorithm AuxiliaryPaths (Algorithm [1]) . 
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Algorithm 1: AuxiliaryPaths(G) 
Input: Graph G = {V, E, ro) 

Output: Set of auxiliary paths for all nodes in G 

II Zi is a 7i-cover of G, < i < ft (ji , k specified in Section [5]) 
// Assume each cluster has a designated leader node 

Q <— 0; // set of auxiliary paths for all nodes in G 
foreach v £ V do 

g(v) <— 0; // auxiliary path for node v from level to s 

x <— v; 

for i = to k — 1 do 

Let X £ Zi+i be a cluster that ^i+i -satisfies v; 
£i+i(v) «— leader of X; 
qi(v) <— shortest path from x to £i+i(v); 
q(v) <— concatenate q(«) and <7i(v); 

return Q; 



Both algorithms use k + 1 covers . . . , Z K , where in Zq every node in V 
is a cluster, and Zi is a 7,-cover of G, for i > 1, where the parameters 7^ and k 
arc defined in Section [5j We refer to the cover Zi as the level i cover of G. We 
assume that each cluster in the covers has a designated leader node. There is a 
unique cluster, containing all nodes in G, and leader node i K at level k. 

Algorithm AuxiliaryPaths computes a auxiliary path q(v) from every node 
v £ V to £ K . The auxiliary paths are built in a bottom- up fashion. A auxiliary 
path from any node v £ V at level is built recursively. In the basis of the 
recursion, we identify a cluster Xi G Zx, which 71-satisfies node v. Let £i(v) 
denote the leader X\. We now compute a shortest path, denoted qo(v), from v 
to tx(v). This forms the first path segment of q(v). Suppose we have computed 
q(v) up to level i, i < k. We now want to extend this path to the next higher 
level i + 1. To compute the path segment from level i to level i + 1, we repeat the 
process of finding a cluster X i+ i 6 Z i+ i that 7^+1 -satisfies node v. Let £ i+ i(v) 
denote the leader X\. We compute the shortest path, denoted qi(v) from £i(v) 
to £i+x(v). We then append this new path segment qt(v) to q(v) to form the 
current extended path q(v). The path building process terminates when the last 
leader reaches level k. 

We are now ready to describe how Algorithm Find Path computes the shortest 
paths between all pair of nodes in G. For pair of nodes u,v £ V, let y be 
the distance between them. Let 7^ be the smallest locality parameter such that 
7i > 2//2. Let X e be the cluster that 7^ satisfies u, and let £' be the respective 
leader of X. Note that by the way that we have chosen 7^, cluster X also ji/2- 
satisfies v. Let q^u) denote the segment of the auxiliary path q(u) from u to 
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Algorithm 2: FindPaths(G) 
Input: Graph G = (V, E, ro). 

Output: Set of paths between all pair of nodes in G 
I / Let covers Zi be as in Algorithm AuxiliaryPaths 

// Let leader li(x) be as defined in Algorithm AuxiliaryPaths; lo(x) <— x 
II Let auxiliary path q(x) be as computed in Algorithm AuxiliaryPaths 
// Let <l(ic) be the auxiliary path segment from x to £i(x); q (x) <— X 

P <— 0; // set of paths for all pairs of nodes in G 
foreach pair u, v 6 V do 

Let 7i be the smallest locality parameter such that 7,: > 2dist(w, v); 

Let X £ Zi be a cluster that 7i-satisfies u (and hence 7i/2-satisfies v); 

Let I be the leader of X {common leader of u, v); 

q' concatenate shortest paths from li~\(u) to £' to li~.\(v)\ 

p(u,v) <— concatenate g i _ 1 (w), <?', and 
_ P ^PUp(u,v); 

return P; 



^i(u). We concatenate q i _ i (u), with a shortest path from ti^i(u) to with a 
shortest path from t' to and This gives the path p(u,v). 

5 Analysis 

Let G = (V, E, to) be a planar graph with n nodes. In this section we use the 
following parameters: 



Consider k + 1 levels of covers Zq, ■ ■ ■ , Z K +\, where in Zq each node in V is a 
cluster, and each Zi, i > 1, is a 7i-cover of G which is obtained from Theorem[TJ 
Thus, each Zi, i > 1, has degree at most f3, stretch at most a, and can be given 
a valid distance^ coloring with \ colors. 

Let A denote an arbitrary set of demands. For any demand d = (s,t) 6 A 
let p(d) = p(s,t) be the path given by Algorithm FindPaths. Suppose that the 
common leader of s and t is t. The path p{d) consists of two path segments: 
the source path segment p(s), from s to I, and the destination path segment p(t) 
from I to t. We denote by Pi(s) the subpath between level i and level i + 1 (we 
call this the level i subpath). 

Let C*{A) denote the cost of optimal paths in A. Let C(A) denote the 
cost of the paths given by our algorithm. We will bound the competitive ra- 
tio C(A)/C*(A). For simplicity, in the approximation analysis, we consider only 



X 



17 



K 



1 + [log 4cr D~\ J /highest cluster level in G 

18 //cover degree bound 

24 // cover stretch bound 

(4cr) l_1 / /locality parameter of level i > 1 cover 

18 //coloring of level i 
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the cost of the source path segments p(si). When we consider the destination 
segments the approximation ratio increases by a factor of 2. 

The cost C(A) can be bounded as a summation of costs from the different 
levels as follows. For any edge e let ip e ^(A) = {pi(s) : ((s, t) £ A)A(e £ Pi{v))} be 
the set of layer-z subpaths that use edge e. Denote by C e ,i(A) = f(\tp e j(A)\) ■ w e 
the cost on the edge e incurred by the level- i subpaths. Since / is subadditive, we 
get C e (A) < J2iZo Ce,i{A). Let C\(A) = J2 e eE C e ,i(A) denote the cost incurred 
by the layer-i subpaths. Since C(A) = J2 eeE ^e(A), we have that: 

C(A)<£ci(^)- (1) 

i=0 

For any cluster X let X(A) denote the set of demands with source in X whose 
paths leave from the leader of X toward the leader of a higher level cluster. 

Lemma 4. For any Z i; 2 < i < n - I, C*(A) > R(i)/\, where R(i) = 
/(|X(A)|). 7l /2. 

Proof. Let Zi(k) to be the set of clusters at level i which receive color k £ [1, %]. 
Consider a cluster X £ Zi(k). Consider a demand (s, t) G X(A). Since X £ Zi(k) 
the common leader of s and t is at a level i + 1 or higher. From the algorithm, 
dist(s,f) > 7j + i/2. Consider the subpaths from X(A) of length up to 7»/2. In 
the best case, these subpaths from X{A) may be combined to produce a path 
with smallest possible total cost f(\X{A)\) ■ ji/2. Any two nodes u £ X(A) and 
v £ X(Y), where X,Y £ Zi(k) and 1^7, have dist(u, v) > 7», since each node 
is 7 i-satisfied in its respective cluster and X and Y receive the same color in 
the distance- 7 i coloring of Z. Therefore, the subpaths of lengths up to 7i /2 from 
the demands X{A) and Y{A) cannot combine. Consequently, C*(A) > R(i,k) 
where R(i,k) = Y^xez^k) f(\ x ( A )\) ' 7i/ 2 - Let i? max = max fee[1;X ] R(i,k). We 
have that C*(A) > i? max . Since R(i) = J2k=i^^) - ^max ■ X- We obtain 
C*(A) > R(i)/x, as needed. 

We also get the following trivial lower bound for the special case where < 
i < 1, which follows directly from the observation that each demand needs to 
form a path with length at least 1. 

Lemma 5. For any Z u < i < 1, C*{A) > J2xez t f(\ x ( A )\)- 
We obtain the following upper bound. 

Lemma 6. For any Z l; < i < k-1, d(A) < Q(i) where Q(i) = Y,xez t f(\ x ( A )\)- 
/3aj l+1 . 

Proof. For any cluster X £ Zi, we can partition the demands X{A) = Y\ U Y2 U 
. . . U Yfe, where Yi j^Yj, i 7^ j, according to leader at level i + 1 that they use, so 
that all demands in Yi use the same leader in Zj+i, and Yi and Yj use a different 
leader of i + 1. Next, we provide a bound on k. 
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Consider any two demands d\ — (si,ii) G X(A) and e?2 = (52,^2) G X(A). 
Let li be the common leader at level i. Since si and s 2 are 7i-satisfied by the 
cluster of li, they are both members of that cluster. Therefore, dist(si,£,) < 
(J7i, and dist(s2, elh) < o^i- Thus, dist(si,sa) < 2aji = 7^+1/2. Suppose that 
demand d\ chooses leader £i + i at level i + 1 with respective cluster Since 
si is at least 7i+i/2-satisfied in Xi+i, s 2 is a member of Xi+i. Since any node 
is a member of at most j3 clusters at level i + 1, it has to be that the number 
of different level i + 1 leaders at level i + 1 that the demands in select is 

bounded by /3. Consequently, k < (i. 

Since / is subadditive and for any demand (s,t) \pi(s)\ < 177^+1, Ci(Yj) < 
/(I^D-^+r. Therefore, C % {X{A)) < Ej=i CiOjO < /(^(^D-^+l Which 
gives: Cj(A) < Lxez, /(l^(^)l) • P°H+ii as needed. 

Lemma 7. For any < i < k - 1, C t {A) < C*(A) ■ 8j3<j 2 x- 

Proof. From Lemma HI for any < i < k — 1, Ci(A) < Q(i). From Lemma [4j 
for any 2 < i < k - 1, C*(A) > R(i)/x- Note that Q(i) = R(i) ■ 2j3axit+illt = 
R(i) ■ 8/3a 2 \- Therefore, C l {A) < C*(A) ■ 8/3a 2 x. For < i < 1, we use the lower 
bound of Lemma O and we obtain < C*(A) • /3aj 2 = C*(A) ■ A/3a 2 . 

We now give the central result of the analysis: 

Theorem 2. The oblivious approximation ratio of the algorithm is O(logn). 

Proof. Since the demand set A is arbitrary, from Lemma [7] and Equation [T] we 
obtain oblivious approximation ratio bounded by 8k/3ct 2 x- When we take into ac- 
count the source path segments together with the destination path segments, the 
approximation ratio bound increases by a factor of 2, and it becomes 16k/3(t 2 x- 
Since, /3, a, Xi are constants and k = O(logn), we obtain approximation ratio 
O(logn). 

6 Conclusions 

We provide a set of paths for the multi-sink buy-at-bulk network design problem 
in planar graphs. Contrary to many related work where the source-destination 
pairs were already given, or when the source-set was given, we assumed the 
obliviousness of the set of source-destination pairs. Moreover, we considered an 
unknown fusion cost function at every edge of the graph. We presented nontrivial 
upper and lower bounds for the cost of the set of paths. We have demonstrated 
that a simple, deterministic, polynomial-time algorithm based on sparse covers 
and shortest-path clustering can provide a set of paths between all pairs of 
nodes in G that can accommodate any set of demands. We have shown that this 
algorithm guarantees 0(log n)-approximation. As part of our future work, we 
are looking into obtaining efficient solutions to other related network topologies, 
such as minor-free graphs. 
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Appendix: Proofs 

Proof of Lemma[J] The algorithm divides the pathp into consecutive disjoint sub- 
paths pi , P2 , ■ ■ ■ , pe each of length 47 (except for the last subpath pi which may 
have shorter length). The algorithm builds a cluster Xi around each subpath 
which consists of the 47-neighborhood of pi. We can show that dist(Xj, > 
7. Suppose otherwise. Then, there are nodes u G Xi, v G which are 7- 

satisfied in their respective clusters and dist(w,w) < 7. Thus, u G Xi+3. Then, 
there is a path of length at most 87 that connects the two paths pi and Pi+3 
which is formed through u. However, this is impossible since the paths are at 
distance at least 87 + 1. Therefore, we can use a palette of at most 3 colors to 
color the clusters, so that each cluster receives color [i mod 3) + 1. 

Proof of Lemma [H Consider a level i > of T. From Lemma [TJ the clusters 
produced in any node w of level i from path p(w) can be colored with 3 colors. 
Consider now two nodes wi and W2 in level i of T. Let X be a cluster from 
p(wi). Any node which is 7-satisfied (with respect to G) in X, cannot have 
in its 7-neighborhood any node in G{w2), since G(w\) and G(w2) are disjoint. 
Therefore, any two nodes which are 7-satisfied in the respective clusters of w\ 
and u>2 have to be at distance more than 7 from each other in G. This implies 
that we can use same palette of 3 colors for each node in the same level i of the 
tree. 

Proof of Lemma [3] From Lemma [21 the clusters of each level of the tree T can 
be colored with 3 colors. From the proof in [5J Lemma 5.4], any node v G G is 
clustered in at most 2 consecutive levels i,i + 1 of T, and does not appear in 
any subsequent level. Any node u which is 7-satisficd in a cluster of level i + 2 
cannot be with distance of 7 or less from v, since v doesn't appear in the level 
i + 2 subgraph of G. Therefore, any node v which is 7-satisfied in a cluster of 
level i must be at distance more than 7 than any node u which is 7-satisfied in 
a cluster of level i + 2. Therefore the clusters formed at level i are at distance 
at least 7 + 1 from clusters formed at level i + 2. Consequently, we can use color 
palette [1,3] for odd levels and color palette [4,6] for even levels, using in total 
6 colors. 



